Observation of the decay B° — > DfAp 
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We report the first observation of the decay B° — > Df Ap with a statistical significance of 6.6 a. 
We measure B(B° -> D+Ap) = (2.9 ± 0.7 ± 0.5 ± 0.4) x 10~ 5 , where the first error is statistical, 
the second is systematic and the third error comes from the uncertainty in B(Df — > 4>ir + ). The 
data used for this analysis was accumulated at the T(45') resonance, using the Belle detector at the 
KEKB asymmetric-energy e + e~ collider. The integrated luminosity of the data sample is 414 fb~ , 
corresponding to 449 x 10 B BB pairs. 

PACS numbers: 13.25.Hw, 13.30.Eg, 14.40.Lb, 14.20.Jn 



In the past few years, new measurements of baryonic 
B meson decays by Belle [1 E RM, @ and CLEO || Q 
have revived experimental [3, [9j, [lO | and theoretical inter- 
est [H E H El EI 13 in such processes. Multi-body 
baryonic decay modes are found to have larger branching 
fractions than two-body modes, and the baryon-pair in- 
variant mass spectrum peaks near threshold in the case 
of multi-body decays [l7|. This feature was conjectured 
in Ref. [3] • Further investigations of the Dalitz plot [n| 
and the angular correlations for events in the threshold 
region [20| offer better understanding of the underlying 
dynamics. 

To date, nothing is known experimentally about 
charmful baryonic B decays with the creation of an 
ss pair. B° mesons can decay to D+Ap through the 
Cabibbo favoured b — > cud process. They can also 
decay to the charge conjugate final state through the 
Cabibbo suppressed b — ► ucd process, opening a new av- 
enue for future CP asymmetry studies. We report here 
the first observation of the decay B° — > DfAp using 
4141b" 1 of data, corresponding to 449 x 10 6 BB pairs, 
collected at the T(4S) resonance with the Belle detec- 
tor at the KEKB asymmetric-energy e + e~ collider [2lj |. 
Since the DfAp final state may get a contribution from 
the D°p — ► DfA final state rescattering, the previously 
observed B° — > D°pp decay [B| could be one of the sources 
for the D+Ap final state. Inclusion of charge conjugate 
states is implicit throughout this paper. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cherenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 



(TOF), and an electromagnetic calorimeter (ECL) com- 
prised of CsI(Tl) crystals located inside a superconduct- 
ing solenoid coil that provides a 1.5 T magnetic field. An 
iron flux-return located outside the coil is instrumented 
to detect K\ mesons and to identify muons (KLM) . The 
detector is described in detail elsewhere [22[ • Two differ- 
ent inner detector configurations were used. For the first 
sample of 152 million BB pairs, a 2.0 cm radius beampipe 
and a 3-layer silicon vertex detector were used; for the 
latter 297 million BB pairs, a 1.5 cm radius beampipe, a 
4-layer silicon detector and a small-cell inner drift cham- 
ber were used We use a GEANT-based Monte Carlo 
(MC) simulation to model the response of the detector 
and determine the efficiency [24I ]. 

Pions, kaons and protons are identified using a like- 
lihood ratio method, which combines information from 
the TOF system and ACC counters with dE/dx mea- 
surements using the CDC [25[ . 

In this analysis we reconstruct Dj candidates by using 
D+ -> K*°K+ and K S K+ decay modes. Candi- 

date A baryons are reconstructed via the A — > p7r~ decay. 

For A hyperons we require an invariant mass within 
±3 MeV/c 2 of the nominal A mass [26|. The distance 
between the A decay vertex position and beam interac- 
tion point (IP) in the r ~ (j) plane, dr(A), is required to 
be greater than 0.5cm. The angle «a> between the A 
momentum vector and the vector pointing from the IP 
to the decay vertex, must satisfy coscva > 0.95. We 
also require dz(A) < 0.5 cm, where dz(A) is the dif- 
ference in the z-coordinates (the z axis is parallel to 
the e + beam) between the 7r and p tracks at vertex 
position. We reconstruct neutral kaons via the decay 
Kg — * 7r + 7r~ and require its invariant mass to be within 
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±10 MeV/c 2 (about 4cr) of the nominal K% mass. We 
also require dz(Kg) < 1cm, dr(Kg) > 0.01cm and 
|cos(a^o)| > 0.95, where dz, dr and a are defined in 
a way similar to the case of the A hyperon. 

We use a mass and vertex constrained fit for Df — > 
K + K ~7r + , and require the 4> invariant mass to be within 
±10MeV/c 2 and K*° invariant mass within ±50MeV/c 2 
of the nominal masses for the — > cf)ir + and £)+ — > 
K*°K + , respectively. Finally, we apply helicity require- 
ments: | cos 0^1 > 0.3 and | cos 67^.0 1 > 0.3 for Df — > 07r 
and Z)+ — > K*°K + , respectively. The helicity angle 
®4>(k*) is defined as the angle between the K(tt) me- 
son momentum and the D s meson momentum in the <f> 
(K*) rest frame. For D s candidates we use a mass win- 
dow that extends ±15MeV/c 2 around the nominal D s 
mass value. We use a large sample of inclusive A and 
D s signals, applying the selections described above, to 
verify that their mass peaks are well described by two 
Gaussians, corresponding to the core and the tail of the 
distribution, where the tail fraction is 35 to 50%. The 
signal mass windows that are used in the analysis corre- 
spond to approximately 4 a for the core and 2 a for the 
tail Gaussian. For the inclusive signals data and MC 
agree. 

To suppress the continuum background (e + e _ — > qq, 
where q = u,d,s,c), we require the ratio of the second 
to zeroth Fox- Wolfram moment [27j to be less than 0.5. 
We also require the cosine of the reconstructed B me- 
son direction with respect to the z-axis in center-of-mass 
(c. m.) frame, | cos6>b|, to be less than 0.8. 

The B candidates are identified by their mass dif- 
ference, AM = M(B) — rriB and their beam-energy 
constrained mass, M bc = x/^bcam ~~ (J2iPi) 2 > where 
-E-beam = \/s/2 is the beam energy and pi are three- 
momenta of the B candidate decay products in the c. m. 
system, M(B) is the reconstructed mass of the B can- 
didate and raj is the world average B meson mass. 
We do not use the widely applied kinematic variable 
AE = Eb — -Ebcam, where Eb is the energy of the re- 
constructed B in the c. m. system, since AE has a 
large correlation with Mb c for signal due to the small 
energy release in the decay under study. By contrast, 
AM and M\, c are uncorrelated (jjl . |29| as confirmed by 
MC. We select B candidates with M bc > 5.2 GeV/c 2 and 
|AM| < 0.2GeV/c 2 . _ 

The M bc and AM distributions for the B° -> D+Ap 
candidates are shown in Fig.[T] (a) and (b), respectively, 
where all three D s modes are combined. We require 
M bc > 5.272 GeV/c 2 (|AM| < 0.025 GeV/c 2 ) for the 
AM (Mbc) projection. We found that after applying all 
the selection requirements, there are no events counted 
repeatedly in the Mb c and AM distributions. The 
hatched histograms in Fig.Q] (a) and (b) show normal- 
ized D s mass sidebands where no peaking structures are 
evident. The superimposed curves are the results of a si- 
multaneous two-dimensional binned maximum likelihood 
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FIG. 1: (a) The M hc and (b) AM distributions for the B° -> 
DfAp candidates (triangles with error bars). The hatched 
histograms show the mass sidebands normalized to the 
signal region. The overlaid curves are fit results (see text for 
details), (c) ThepA and (d) DfA invariant mass distributions 
in the B-signal region (open histogram) and in the B-sideband 
(hatched histogram). 



fit (with common branching fraction as a constraint) to 
the three AM versus Mb c distributions (for the three D s 
channels). 

To describe the signal we use Gaussians with means 
and widths fixed to the values obtained from MC. The 
backgrounds in Mb c and AM are parameterized by a 
first-order polynomial and an ARGUS function [301 ] . re- 
spectively. The fit gives a statistical significance of 6. Go- 
fer the signal, where the statistical significance is defined 
as yf— 2 ln(Lo/^max)j where Lq and £ ma x are the likeli- 
hoods with the signal fixed at zero and at fitted value, 
respectively. The region AM < — 0.08GeV/c 2 is ex- 
cluded from the fit to avoid possible contributions from 
the B° -> D*+Ap, D*+ -> D+7 and B~ -> D+Apn~ 
decays, where the soft 7(71"^) is undetected. The choice 
of the fitting range is taken into account in the system- 
atic error. The results of the fit applied for the three 
modes separately are shown in TableQ] 

We select events in the B-signal region of |AM| < 
0.025 GeV/c 2 and M bc > 5.272 GeV/c 2 for the three D+ 
modes and examine the two-baryon invariant mass dis- 
tribution (Fig.[T](c)). We see apparent threshold peaking 
behavior of this distribution, while the -B-sideband [31( 
distribution is smooth. Such a threshold peaking behav- 
ior seems intrinsic to all multi-body baryonic B decays. 
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TABLE I: Summary of the fit results, efficiencies, statistical significances and branching fractions obtained from the 2D 
AM - M bc fit. 



Decay Mode 


Yield 


Efficiency (10~ 4 ) 


Significance 


B (10~ b ) 


B° -» DfAp, Df -> (f>-K + 


6.5 ±2.6 


4.90 


4.7(7 


3.0 ± 1.2 


B° -> DfAp, Df -> K*°K + 


4.0 ±2.5 


4.31 


2.3cr 


2.1 ± 1.3 


B° -» DfAp, Df -> K° S K+ 


7.9 ±3.1 


4.83 


4.2cr 


3.6 ± 1.4 


B u — > DfAp, simultaneous fit 






6.6(7 


2.9 ±0.7 
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FIG. 2: Cross-checks of the signal, (a) and (b) are the Mbc 
and AM distributions for the on-resonance data with the in- 
verted requirement on the normalized Fox- Wolfram moment 
i?2 (see text for details), (c) and (d) are the Mb c and AM 
distributions for the continuum data, (e) and (f ) are the Mb c 
and AAf distributions for the primary vertex protons with 
inverted identification. No peaking structure in the signal 
region is present in any of the cross-check analyses. 



The invariant mass distribution for Df and A and for 
corresponding .B-sideband is represented as well (Fig.[T] 
(d)). Some peaking behavior is also seen in this distri- 
bution, which could arise from some new excited charm 
baryon. Firm conclusions, however, cannot yet be drawn 
on either peaks because of limited statistics. 

As a cross-check we analyze the on-resonance data with 



the inverted requirement on the normalized Fox- Wolfram 
moment i?2 > 0.5 (Fig. [2] (a), (b)) and off-resonance data 
sample (Fig. [21(c), (d)) and find no candidates in the B- 
signal region. The distributions for the primary vertex 
protons with an inverted particle identification require- 
ment [sH do not peak in the signal region (Fig. [21(e), (f)), 
demonstrating that the selected B° — * DfAp candidates 
contain real protons. 

Tablcfl] summarizes the results of the fits 
construction efficiencies including the B(Df 
B{(j> -> K+K-), B{Df -> K* a K+) 1 B(K*° > 



B(Dj 



K° S K+), B(K° S -> 7T+7T-), B(A 



the re- 

- 07T+), 
K-TT+), 



branching fractions, statistical significance of the signals 
and extracted branching fractions. Here we assume equal 
fractions of charged and neutral B mesons produced in 
T(4S") decays. 

The major sources of systematic error are the uncer- 
tainties in the tracking efficiency 6% (1% per track), 12% 
in the charged particle identification efficiency (1% for 
pion, 2% for kaon, 3% for proton), 5% for A finding, 3% 
for efficiency estimation due to MC statistics and 5% for 
the error due to choice of the fitting procedure. These 
contributions are combined in quadrature resulting in a 
total systematic error of 16%. We also take into account 
a third error due to the uncertainty in B(Df — ► <j)-K + ) 
that is 14%. 

In summary, we report the first observation of the de- 
cay B° — > DfAp with a branching fraction of (2.9 ± 
0.7±0.5±0.4) x 10 -5 , where the first error is statistical, 
the second is systematic, and the third arises from un- 
certainty in the branching fraction of Df — > cf>ir + . The 
statistical significance is 6.6 a. This charmful decay can 
occur via the creation of an ss pair or from the more 
copious B — > DNW modes with (DN) + — * DfA rescat- 
tering in the final state [B|, 0, HI]- In the future, this 
decay mode can be used for CP asymmetry studies. 

We thank the KEKB group for excellent operation of 
the accelerator, the KEK cryogenics group for efficient 
solenoid operations, and the KEK computer group and 
the Nil for valuable computing and Super-SINET net- 
work support. We acknowledge support from MEXT and 
JSPS (Japan); ARC and DEST (Australia); NSFC and 
KIP of CAS (China); DST (India); MOEHRD, KOSEF 
and KRF (Korea); KBN (Poland); MES and RFAAE 
(Russia); ARRS (Slovenia); SNSF (Switzerland); NSC 



■5 



and MOE (Taiwan); and DOE (USA). 



[1] Belle Collaboration, N. Gabyshev et al, Phys. Rev. 

Lett. 90, 121802 (2003). 
[2] Belle Collaboration, M.-Z. Wang et al, Phys. Rev. 

Lett. 90, 201802 (2003). 
[3] Belle Collaboration, K. Abe et al, Phys. Rev. Lett. 88, 

181803 (2002). 

[4] Belle Collaboration, N. Gabyshev et al, Phys. Rev. D 

66, 091102 (2002). 
[5] Belle Collaboration, K. Abe et al, Phys. Rev. Lett. 89, 

151802 (2002). 

[6] CLEO Collaboration, S. A. Dytman et al, Phys. Rev. D 

66, 091101 (2002). 
[7] CLEO Collaboration, S. Anderson et al., Phys. Rev. 

Lett. 86, 2732 (2001). 
[8] Belle Collaboration, M.-Z. Wang et al, Phys. Rev. 

Lett. 92, 131801 (2004). 
[9] Belle Collaboration, K. S. Park, H. Kichimi et al., Phys. 

Rev. D 75, 011101 (2007). 
[10] Belle Collaboration, M.-Z. Wang, Y.-J. Lee et al, Belle 

preprint 2007-19, arXiv: 0704.2672[hep-ex]. 
[11] C.-K. Chua, W.-S. Hou, S.-Y. Tsai, Phys. Rev. D 65, 

034003 (2002). 

[12] C.-K. Chua, W.-S. Hou, S.-Y. Tsai, Phys. Lett. B 528, 
233 (2002). 

[13] J. L. Rosner, Phys. Rev. D 68, 014004 (2003). 

[14] B. Kerbikov, A. Stavinsky and V. Fedotov, Phys. Rev. C 

69, 055205 (2004). 
[15] J. Haidenbauer, Ulf-G. Meissner and A. Sibirtsev, Phys. 

Rev. D 74, 017501 (2006). 
[16] D. R. Entem and F. Fernandez, Phys. Rev. D 75, 014004 

(2007). 

[17] H. Kichimi, Nucl. Phys. B Proc. Suppl. 142, 197 (2005). 
[18] W.-S. Hou, A. Soni, Phys. Rev. Lett. 86, 4247 (2001). 
[19] BABAR Collaboration, B. Aubert et al, Phys. Rev. D 



72, 051101 (2005). 
[20] Belle Collaboration, M.-Z. Wang et al, Phys. Lett. B 

617, 141 (2005). 
[21] S. Kurokawa and E. Kikutani, Nucl. Instr. Methods Phys. 

Res., Sect. A 499, 1 (2003), and other papers included 

in this volume. 

[22] Belle Collaboration, A. Abashian et al, Nucl. Instr. 
Methods Phys. Res., Sect. A 479, 117 (2002). 

[23] Z. Natkaniec et al. (Belle SVD2 Group), Nucl. Instr. and 
Meth. A 560, 1 (2006). 

[24] R. Brun et al, GEANT 3.21, CERN DD/EE/84-1, 1984. 

[25] Charged kaons are required to satisfy C{K)/{C(K) + 
£(7r)) > 0.6. Charged pions are required to satisfy 
C(n) / (C(K)+£.(tt)) > 0.1. Protons are required to satisfy 
C{p)/{C{K)+C{p)) > 0.6 and£(p)/(£(7r)+£(p)) > 0.6. 
Here C(K/n/p) is the particle identification likelihood for 
the K/tt/p hypotheses. The above requirements have ef- 
ficiencies of more than 95% for pions, kaons and protons, 
respectively, from B° — ► D+Ap decays. The probability 
for each particle species to be misidentified as one of the 
other two is less than 5%. 

[26] W.-M. Yao et al, (Particle Data Group) J. Phys. G 33, 
1 (2006). 

[27] G. C. Fox and S. Wolfram Phys. Rev. Lett. 41, 1581 
(1978). 

[28] Belle Collaboration, S. L. Zang et al, Phys. Rev. D 69, 
017101 (2004). 

[29] Belle Collaboration, N. Gabyshev et al, Phys. Rev. 
Lett. 97, 202003 (2006). 

[30] ARGUS Collaboration, H. Albrecht et al, Phys. Lett. B 
241, 278 (1990). 

[31] The £?-sideband is defined as a region of Mb c > 
5.2 GeV/c 2 and -0.08 GeV/c 2 < AM < 0.12 GeV/c 2 ex- 
cluding the B-signal region. 

[32] Protons are required to satisfy inverted likelihood ratio 
requirement £{p)/(£(K) + C(p)) < 0.6. 

[33] BABAR Collaboration, B. Aubert et al, Phys. Rev. D 
74, 051101 (2006). 



